Spin-echo dynamics of a heavy hole in a quantum dot 
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We develop a theory for the spin-echo dynamics of a heavy hole in a quantum dot, accounting 
for both hyperfine- and electric-field-induced fluctuations. We show that a moderate applied mag- 
netic field can drive this system to a motional-averaging regime, making the hyperfine interaction 
ineffective as a decoherence source. Furthermore, we show that decay of the spin-echo envelope is 
highly sensitive to the geometry. In particular, we find a specific choice of initialization and 7r-pulse 
axes which can be used to study intrinsic hyperfine-induced hole-spin dynamics, even in systems 
with substantial electric-field-induced dephasing. These results point the way to designed hole-spin 
qubits as a robust and long-lived alternative to electron spins. 
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Electron spins in solid-state systems provide a versa- 
tile and potentially scalable platform for quantum infor- 
mation processing [l], Q. This versatility often comes 
at the expense of complex environmental interactions, 
which can destroy quantum states through decoherence. 
Many theoretical and experimental studies have now es- 
tablished that the coherence times of electron spins in 
quantum dots 0- (1 , bound to donor impurities 0, Q , and 
at defect centers 9( are typically limited by the stron g hy - 
pcrfine interaction with surrounding nuclear spins [2tllfj|. 
Heavy-hole spin states in III-V semiconductor quantum 
dots have emerged as a platform that could mitigate 
the negative effects of the hyperfine interaction. Due 
to the p-like nature of the valence band in III-V ma- 
terials, the contact interaction vanishes for hole spins, 
leaving only the weaker anisotropic hyperfine coupling 
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12j . Moreover, the anisotropy of this interaction in 
two-dimensional systems should allow for substantially 
longer dephasing times in a magnetic field applied trans- 
verse to the quantum-dot growth direction 10, 11 1 . 



Recent experiments have measured hyperfine coupling 
constants for holes fl3T - fl5^ . as well as spin-relaxation (T\) 
[HI and free- induction decay times, T 2 *, through indirect 
(frequency-domain) [I?} and direct (time-domain) stud- 



ios 
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Coherent optical control has now been demon- 
strated for hole spins in single [3, [lj| and double quan- 
tum dots This technique has been used to imple- 



ment a Hahn spin-echo sequence [18l ] giving an associated 
spin-echo decay time, T 2 ~ 1 /is. The T 2 value reported 
in Ref. [l8| has been attributed to device-dependent 
electric-field fluctuations, rather than the intrinsic hyper- 
fine interaction. Motivated by these recent experiments, 
here we present a theoretical study of heavy-hole spin- 
echo dynamics with an emphasis on identifying the opti- 
mal conditions for extending coherence times. In partic- 
ular, we show that dephasing due to electric-field fluctua- 
tions, as proposed in Ref. |18(, is dramatically suppressed 
in an alternate geometry considered here. Moreover, in 
contrast with the case of electron spins, we find that hole 
spins can enter a motional-averaging regime in a moder- 
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FIG. 1. (Color online) (a) Quantum-dot geometry, with nu- 
clear field h z and magnetic field B = Bx. For unstrained and 



flat quantum dots (d < L), 7// 



_/zb ^ and h x 







[111. |21|. |22||. (b) Hahn echo sequence. Two 7r-rotations (Un 
and Ul, taken here about the ic-axis) are applied at t = t and 
2t, to refocus the HH spin. 



ate magnetic field. In this regime, coherence is no longer 
limited by the hyperfine interaction, solidifying the po- 
tential for long-lived hole-spin qubits. 

We consider a heavy-hole (HH) spin interacting with 
nuclear spins in a flat quantum dot with weak strain (see 
Fig. [T]). The HH spin is then described with the following 
Hamiltonian 23[ (setting h = 1), 



H = H P 



h z S z 



- 1H BS X - 



BI, 



k ■■ 



(1) 



where S = cr/2 is a pseudospin-1/2 operator in the two- 
dimensional (J z = ±3/2) HH subspace and Ik the nu- 
clear spin at site k. Hz gives the hole- and nuclear- 
Zeeman interactions for an in-plane magnetic field B = 
Bx [see Fig. Ha)]. li k i s the gyromagnetic ratio of iso- 
tope ik at site k having total spin Ii k . The hole gyromag- 
netic ratio is = g±^B, with g± the in-plane g- factor 
for a dot with growth axis along [001] and [Ib the Bohr 
magneton. The hyperfine interaction [10l.lll| is expressed 
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in terms of the Overhauser operator, h = V^ fe A^I^. The 
coupling constants, Ak, are given by A k = A lk v \ip(rk)\ 2 , 
with A 1 the hyperfine constant for isotope i, Vq the vol- 
ume occupied by a single nuclear spin, and ip( r k) the 
HH envelope wavefunction. When the isotopes are dis- 
tributed uniformly across the dot, we define the average 
A = ^2 k Ak ~ £\ viA % , with vi the isotopic abundance. 
In this case, and for a Gaussian envelope function in two 
dimensions, A k ~ [A/N)e~ k / N [3] with TV = 10 4 - fO 6 
a typical number of nuclear spins within a quantum-dot 
Bohr radius. The ratio of \A\ to the strength of the hy- 
perfine coupling of electrons, has been estimated 
theoretically in GaAs and confirmed experimentally 
[II 01] in InGaAs and InP/GalnP to be {A/A^l ~ 0.1. 
For simplicity, we will evaluate numerical estimates with 
a single averag ed value \A l \ ~ \A\ ~ 13 peV [l^, EH and 
appropriate for Ino.5Gao.5As. 
Spin echo. Under the action of Eq. ([T]), spin dephas- 
ing results from fluctuations in h z . Provided these fluc- 
tuations remain static on the timescale of decay of the 
hole spin, this source of decay can be removed via a Halm 
echo [see Fig.QJb)]. The process is better analyzed in the 
interaction picture with respect to Hz, 



H(t) = h z (t)S z (t), 



where, for any O, 0(t) 



Oe 



iH z t 



(2) 
In par- 



ticular, h z (t) = *£ k A k [q cob (j ik Bt) - 7^sin(7 ife St)] 
and S z (t) = [S z cos{j H Bt) - S y sin{j H Bt)}. The time- 
evolution operator after a time 2r is then given by: 



U(2t) = Te~ i ^ dtfl ^ t) . 
Here, the modified echo Hamiltonian, 



H e (t) 



H(t) 

a a H(t)a a 



< t < T, 

t < t < 2r, 



(3) 



(4) 



takes into account 7r Q -pulses (7r-rotations about a = 
x,y,z). As seen in Eq. (|4|), ^-pulses (but in general 
not TTy , ir z ) have the beneficial effect of inverting the 
sign of the Hamiltonian, Hit) — > —H(t), in the interval 
t < t < 2r. Provided H(t) is approximately static over 
the interval < t < 2t, this will induce time-reversed 
dynamics for r < t < 2t, refocusing decay at the time 
2t. For this reason, unless otherwise specified, we will 
focus in the following discussion on a geometry with the 
magnetic field along x and ^-pulses. We will contrast 
this analysis later with an alternate geometry relevant to 
recent experiments. 

Vanishing g±. We first consider the limit = 
.9j_Mb = in Eq. (fTJ). The dynamics we find in this limit 
will be a good description whenever jh < 7*j correspond- 
ing to q 1 < 10~ 3 (g± < 5 x 10~ 3 has been reported in 2D 
wells |25j ) . This limit considerably simplifies the theoret- 
ical analysis and allows for an exact solution: H becomes 
block diagonal in the eigenbasis of S z and, in each block, 
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FIG. 2. (Color online) Exact analytical spin-echo decay with 
B = 50 mT (solid gray), B = 2 T (dash-dotted), and 73 = 10 T 
(solid black), corresponding to the regimes B < A/^iyN), 
B ~ A/(7,v / TV), and B > A/(<jiVN), respectively. Mark- 
ers show the approximate form, Eq. (Q). We have chosen 
(S x (0)) =1/2, TV = 10 4 , 7 h = g ± p B = 0, and 7i from Table 
1 of Ref. fTfll. 



the eigenstates are obtained after rotating 7^-eigenstates 
by an angle 9^ = ± arctan(27 ifc i3/^4fe) about y. Repre- 
sentative results of the exact evolution of (S x (2t)) are 
shown in Fig. [5J The spin-echo signal has a remarkable 
dependence on the magnetic field: there is a clear tran- 
sition from a low-field regime, where the decay time de- 
creases with increasing B, to a high-field regime, where 
there is no decay, only modulations of the echo envelope. 

To give physical insight, we have developed an analyt- 
ical approximation scheme based on the Magnus expan- 
sion. The Magnus expansion is an averagc-Hamiltonian 
theory typic ally applied to periodic and rapidly oscillat- 
ing systems [26| . This scheme is suggested by the oscillat- 
ing terms in Eq. @ , and will allow us to analyze the more 
general problem with 7^ 0. In the Magnus expansion, 
we assume the evolution operator, Eq. ([3]), can be writ- 



ten as U(2t) 



-%Hm[2t) 



,H^>(2r)_ Thc 



;th 



order term, ff W (£) ; is found using standard methods [2a ]. 
Each higher-order term in the Magnus expansion con- 
tains one additional integral over time. Oscillating terms 
are therefore suppressed by a factor of order ||if ||/u;, with 
oj the typical oscillation frequency. The leading-order 
term is if (0) (i) = H(t)t, where H(t) is simply the aver- 
age of H e (t) over an interval t. The spin components S a 
(a = x, y, z) are then given by: 

(S a (2r)) = (uH2T)S a (2T)U(2T)) = ( e S L »^ S a (2r)) , 

(5) 

where L M {t) is defined by L m (2t)0 = [H m (2t), O] and 
(O) = Tr{Op}. The initial state p = ps ® pi is assumed 
to describe a product of the hole-spin (ps) and nuclear- 
spin (pi) density matrices, where the nuclear spins are 
in an infinite-temperature thermal state. For TV ^ 1 un- 
corrected nuclear spins, the central-limit theorem gives 
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nearly Gaussian fluctuations, resulting in 

(e""™S a ) ~ (exp{ - i(ii / (2r)) / }5 Q 



(6) 



where we define (L 2 M {t))jO s = Tv I {(L 2 M (t)O s )pi} and 
(0) s = Tv s {Op s }. 

At high B, rapid oscillations in H(t) allow us to keep 
only the leading term: L m (2t)G s ~ Z,( o )(2t)0,s = 
[H^(2t), O s ]. Setting 7^ = 0, as appropriate for Fig.0 
and with the help of J2'k A k - ViA 2 /(2N) (where the 
prime restricts the sum to nuclei of isotopic species i), 
we obtain: 



(W) 



cxp 



E 



AvjA 2 !^ + 1) 
3JV( 7< B) a 



(7) 

As seen in Fig. [2] Eq. (J7J (markers) reproduces the ex- 
act dynamics very well. The precise conditions for the 
validity of the Magnus expansion will be given below. 

The simple form of Eq. Q enables us to under- 
stand why the behavior of (S x (2t)) changes as B is in- 
creased. For B <C A/{^/i^fN) (gray solid line in Fig. [2}, 
a short-time expansion of Eq. (J7J gives (S x (2t)) ~ 
(5,(0)) (1 - (r/r ) 4 ) ~ (S x (0))e-^\ with 



1 



E 



^( 7 ,A) 2 ^(Jj + l) 
4iV 3 



-1/4 



(8) 



Surprisingly, when B is increased, r D decreases. This be- 
havior is opposite to the situation for electron spins, in 
which the echo decay time increases for increasing B [27j. 
This decrease is due to rapid fluctuations in h z from nu- 
clear spins precessing at frequencies ~ JiB. The Hahn 
echo can no longer refocus these dynamical fluctuations 
at finite B, although Eq. ([7} does predict partial re- 
currences (dash-dotted line in Fig. [5]) due to the finite 
number of discrete precession frequencies ~ jiB. 

In contrast, at large magnetic field, B > A/^y/N), 
the system enters a motional-averaging regime in which 
the decay of (S x (2t)) is bounded by - (A/ B^y/N) 2 , 
giving rise to beating (black solid curve in Fig. [2]). This 
beating has the same physical origin as electron-spin- 
echo envelope modulation (ESEEM) [HI, although the 
extreme anisotropy of the hole hyperfine interaction al- 
lows uniquely for its complete suppression. Fig. [3] shows 
the 1/e decay time, t q , as B is increased, leading to a dis- 
continuity when jiBy/N J A > 1, at which point (S x (2t)) 
always remains close to its initial value [see Eq. 0]. 

Finite g±. Although there are definite advantages to 
making flat unstrained dots leading to g± ~ and h x , v — 
0, current experiments are performed on hole systems 
with finite (albeit small) g± [3, [2^] . For this general 
case, with 7^ = g±(iB =/= 0, we have no closed-form exact 
solution for the dynamics, but our analysis can still be 
applied for a certain range of r, B. 




FIG. 3. (Color online) Decay time, r , vs. B. Here, 
7 = £\ Vi^fi and parameters are as given in the caption of 
Fig. [2] Insets: typical (S x (2t)) in each of the three regions: 
^fBVN/A < 1, yBy/N/A ~ 1, and ^jByfN/A > 1. 



We neglect subleading oscillating terms in the Magnus 
expansion when ||(# (0) ) 2 || > \\H^H^\\,\\ (H^) 2 \\. 
More specifically, if the relevant fast oscillation frequency 
is uj ~ 7i-Bj each precessing nuclear spin experiences 
a typical hyperfine field <5w rms ~ A/N from the hole. 
Otherwise, if the fast frequency is uj ~ JhB, the hy- 
perfine field acting on the precessing hole is of order 
Su) Tms ~ A/y/N, averaging over the nuclear configura- 
tions. As a consequence, the parameter 5uj Itns /uj < 1 
controls the expansion with Soj Trns = A/y/N (A/N) and 
uj = jhB (jiB) for j H S> 7, (7/r -C ji). In addition to 
bounded oscillating terms, the Magnus expansion gener- 
ates terms that grow with r. These terms approach ~ 1 
at t ~ T max , beyond which a finite-order Magnus expan- 
sion fails. For B — > 00, we find that progressively higher- 
order terms become relevant at a progressively shorter 
time, saturating at r max ~ uj/Sui 2 ms . However, for the 
special case jh -C 7« there is an intermediate regime, 
5uj Tms <C JiB -C y/NSuirms, in which low-order sub- 
leading terms dominate, giving r max ~ (l/A)(uj/Suj rms ) 3 . 
There is no such intermediate regime when 7^ 3> 7i- 
The upshot of this analysis is that uj > 5cj Tms (i.e. 
jiB > A/N) can be easily satisfied for "C ji, even 
outside of the motional-averaging regime (defined by 
7i-B > A/y/N), giving the good agreement shown in Fig. 
[2] In contrast, when jh S> 7i, the Magnus expansion 
is only strictly valid in the motional-averaging regime, 
jhB > A/y/N (Fig. Nevertheless, in this regime the 
Magnus expansion will provide an accurate description 
whenever 8uj Tms /u) < 1 and for r < r max . 

In practice, the value g± = 0.04 measured in (29| sug- 
gests that jh 3> 7i in many current experiments. For 
gj_ = 0.04, the condition B > A/^hVN) is already 
satisfied above a rather small value, B > 100 mT for 
N ~ 10 4 . In Fig. U we plot representative curves in this 
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FIG. 4. (Color online) Main panel: spin-echo envelope from 
a leading-order Magnus expansion when 3> 7» in the 
motional-averaging regime, jhB > A/yN. Inset: Time- 
averaged (S x (2t)). We have taken g± = 0.04 (S x (0)} = 
1/2, and N = 10 4 . These results are valid at least for 
2r < 2r max ~ 70ns (140ns) for B = 400 mT (B = 800 mT). 



motional-averaging regime, displaying the same features 
discussed for g± = 0. Additionally, fast oscillations at 
the hole Zeeman frequency, JhB, induce beating in the 
echo envelope function, (S x (2t)}, which is not present for 
9x = 0. 

Decay anisotropy. While the results discussed so far 
are specific to Tr^-pulses, other schemes are possible. Due 
to the extreme anisotropy of the hole-spin hyperfine cou- 
pling, the spin-echo decay is also highly anisotropic, de- 
pending on both the initialization and 7r-pulse axes. If 
the hole spin is initialized along a generic in-plane di- 
rection h = n x x + n y y, and 7r-rotations are performed 
about that same axis, we find that (S(2r) • h) / (S(0) ■ n) 
is independent of n when g± = 0. This result is to be 
expected since any in-plane component of the hole spin 
experiences the same effective field, h z (t), along the z- 
axis. On the other hand, for g± ^ 0, rotational sym- 
metry about the z-axis is broken, resulting in a strong 
in-plane anisotropy. For the parameters of Fig. 01 but 
with initialization along y and ^-pulses, we obtain that 
(S v (2t)) is dominated by the hole Larmor precession 
about the x-axis and approaches the simple sinusoidal 
function (S v (2t)) ~ (S y (0)) cos(2-fn Bt) in the motional- 
averaging regime, ^hB > A/yfN . 

Additional dephasing mechanisms other than the nu- 
clear bath can also have a strong influence on the pre- 
cession about i, introducing other sources of anisotropy. 
In particular, the decay of (S v (2t)) was measured in 
18[ with 7r z -pulses used for the Hahn echo. The re- 
sulting decay was found to be approximately exponen- 
tial, (S y (2r)) ~ (S y (0))e- 2r ^ T ' 2 , with a B-independent 
T-2 ~ 1 /is. This behavior was attributed to spectral dif- 
fusion induced by electric-field noise, which we model 
here by setting ~f H B ->■ ~f H B + 5oj(t) in Eq. §Q 
The observed exponential decay is consistent with Gaus- 



sian white noise of the form (5w(t)6uj(t'))su = ijrS(t — t') 
(and (5u}(t))su = 0), where (. . .}$u indicates averaging 
with respect to realizations of 6ui(t). We have included 
this additional dephasing mechanism in the evaluation of 
Eq. ([5]) for the 7r K -pulse echo sequence examined previ- 
ously and obtained a power-law decay at r 3> 



(W) 
(SAO)) 



(cxp 



- E m 



i 



> 5u 



1 + t/t d ' 



(9) 

where f y (r) = J T sin 0(i)sgn(r - t)dt, / 2 (r) = 
^ T cos 4>(t)dt, 4>(t) =j H Bt + f* 6w(t')dt' , and 



Td 



l + { lH BT 2 f 
1(hl) x T 2 



(10) 



This decay time scale is exceedingly long (td ~ 20 s) 
for the experimental value jhB ~ 2 x 10 11 s" 1 and us- 



ing (h 2 z 
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is , 



which demonstrates the negligi- 



ble effect of spectral diffusion on the previous discussion 
(e.g., Figs. [2J [31 and SI). For simplicity, we have derived 
Eqs. ([9]) and (fTO)) with static nuclear- field fluctuations 
(h 2 ) i. This corresponds to a worst-case scenario for the 
present model. At the high magnetic field of Ref. [l8| 
(B ~ 8 T), motional averaging would likely inhibit decay 
even further. 

Conclusion. We have calculated the spin-echo dy- 
namics of a single heavy-hole spin in a flat unstrained 
quantum dot. The relevant dynamics are highly 
anisotropic in the spin components and 7r-rotation axes. 
When jh ^Ti, we predict an initial decrease of the co- 
herence time with increasing B, followed by a complete 
refocusing of the HH-spin signal and motional averaging 
when B >B C (B c ~ A/^VN ~ 3T for N = 10 4 ). The 
motional-averaging regime is also realized when 3> 7; , 
relevant to current experiments. In this regime, de- 
cay due to the hyperfine coupling can only occur for 
T ^ Tmax oc B, and can therefore bc completely sup- 
pressed. We have further shown that device-dependent 
electric-field noise becomes negligible for a specific geom- 
etry, allowing for a measurement of the limiting intrinsic 
decoherence due to nuclear spins. 
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